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The thermal nucleation of quark matter bubbles inside neutron stars is examined 
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for various temperatures which the star may realistically encounter during its lifetime. 
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PsJ ' It is found that for a bag constant less than a critical value, a very large part of the 
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^^ \ star will be converted into the quark phase within a fraction of a second. Depending 

,—1. ' on the equation of state for neutron star matter and strange quark matter, all or some 

of the outer parts of the star may subsequently be converted by a slower burning or 
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I ' a detonation. 
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I. INTRODUCTION 

If pulsars or the central parts of these can be made of quark matter rather than neutrons 
|jl|], does this then apply to all or just some of them, and when and how does the phase 
transformation take place? 

According to some investigations 0, the transformation occurs during the supernova 
explosion. In this scenario, the released binding energy is what makes the supernova succeed 
in the first place, supplying the final "push" which seems to lack in most of the computer 
simulations of the events. 

Another model for strange star formation (in the context of absolutely stable strange 
quark matter) was introduced by Baym et al. 0, describing the transformation as a slow 
burning (combustion) rather than a violent event connected with a supernova detonation. 

Regardless of the way in which the transformation occurs, an initial seed of quark matter 
is needed to start it. Alcock et al. [§[ suggested a variety of possibilities ranging from pressure 
induced conversion via two flavor quark matter to collision with either highly energetic 
neutrinos or smaller lumps of strange quark matter. However, they did not provide a rate 
for conversion of neutron stars and thus left it as an open question, whether every compact 
object is a strange star, or whether they are rare objects, even if quark matter formation 
is energetically favorable. It has also been suggested [Q that strange matter seeds (in the 
case of quark matter stability) from strange star collisions or of cosmological origin would 
trigger the transformation of all neutron stars, in which case the thermal nucleation would 
be of relevance to the case of unstable quark matter only P|. 

Other possibilities are that shock waves in the supernova trigger the conversion; a seed 
could be produced by non-thermal quantum fluctuations, or a phase transition could be 
started around impurities. We are not able to estimate the probability of either method, 
but would expect at least quantum fluctuations to be less likely than the thermal nucleation 
process discussed below. 

An estimate for quark matter formed via thermally induced fluctuations was given by 



Horvath et al. 0], using typical numbers for various physical quantities. It was found that 
all neutron stars are converted into strange stars (assuming stable strange quark matter) if 
the temperature at some time during the stars lifetime has exceeded 2-3 MeV. 

In the following, we choose an approach similar to the one in Ref. [0, but with an extra 
term in the expression for the surface energy of quark matter, and with the two phases 
treated in a more self-consistent way. Unlike most of the approaches mentioned above 
|ll|J^J^,|,|],0 , we will be considering the formation of both strange stars (for absolutely stable 
strange matter) and hybrid stars, where strange matter is formed only in the central regions 
due to the high pressure. First, we will deal with some general aspects of nucleation (Sec. 
IID. For pedagogical purposes. Sec. |ITl| treats the problem using a simplified model with the 



hadron phase being a free degenerate neutron gas and the quark phase a bag model with 



only u and d quarks. A more detailed model for the neutron star is presented in Sec. |V 
followed by some concluding remarks (Sec. 0). 



II. BUBBLE FORMATION 

The free energy involved in formation of a spherical quark bubble of radius R is given 

by 

Atx 
F = -—R^AP + 47T(tR^ + 87f-fR + NgAfi, (1) 

where AP = Pq — P^ is the pressure difference, a = aq + ah the surface tension, 7 = 7g — 7/i 
the curvature energy density, and A/i = fiq — fih the difference in chemical potential. A''^ 
is the total baryon number in the quark bubble. The indices h and q denote the hadron 
and quark phase respectively. (We set T = in the thermodynamical expressions; since 
T <^ fi throughout, this only leads to minor errors). Defining C = C{^q) = AP — UqA/j, 
and b = 2'~fC/a'^ the free energy has a maximum at the critical radius 

rc = ^{l + Vl + b), (2) 

and the corresponding free energy 



W, ^ Fir,) = -^ [2 + 2(1 + bf' + 36] (3) 

is the work required to form a bubble of this radius which is the smallest bubble capable of 
growing. 

It is a standard assumption in the theory of bubble nucleation in first order phase tran- 
sitions that bubbles form at this particular radius at a rate given by |^ 

n^T^exp{~WjT). (4) 

It is possible to show that W, has a minimum as a function of fig at fig = /ih- This gives 
a maximum in the rate for bubble formation, and because of the exponential in Eq. (§) one 
may safely assume that nucleation happens in chemical equilibrium. Thus, C reduces to 
C = AP. 

Throughout this paper, we will consider the strange quark to be massless, in which case 
chemical equilibrium would give /i^ = /^d = /is = ^fJ'h- However, equilibrium is established 
only on a weak interaction time scale, whereas the formation of bubbles is governed by the 
strong interaction and is many orders of magnitude faster. So, although we have chemical 
equilibrium between the two phases this is not so between the three quark flavors. Instead, 
flavor must be conserved during the phase transition. 

A consequence of having massless quarks is that cr^ = [^], and since dh is negligible 
compared to 7, Eqs. (0) and (^) reduce to 
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^ ^ AP 
and 



(5) 



»-. = — 1?^. (6) 

3 V AP ^ ^ 



III. PURE NEUTRON GAS 

Before considering a more realistic equation of state it is instructive to study the nucle- 
ation of a pure neutron gas into quarks. The quark bubbles formed consist of u and d quarks 



in the ratio 1:2; only later weak interactions may change the composition to an energeti- 
cally more favorable state. Thus quark chemical potentials are related by /i^ = 2^/^/i„, and 
fi'Ti = f^u + '^l^d = (1 + 2'^/^)/iu, assuming chemical equilibrium across the phase boundary. 
The pressure difference is given by 

AP = P„, - P„ = ^^^i^ - P - P„ (7) 



and the curvature energy coefficient [^ 



7=^. (8) 

For the question in hand we choose the simplest possible equation of state for the neutron 
gas, namely that of a zero temperature, nonrelativistic degenerate Fermi-gas, where 



Ibir'^rrir, 



(9) 



and the baryon density 



ns ^ M^!^ (10, 



A necessary condition for bubble nucleation is that AP > 0. This leads to an upper limit 
on the bag constant, Bmax, from Eq. (|^) as illustrated in Fig. |l| (the corresponding limit for 
the Bethe- Johnson equation of state is shown for comparison; it is seen to be very similar). 

Also shown in Fig. |l]is the limit on the bag constant below which bubble nucleation takes 
place at rates exceeding 1 km^'^Gyr^^ and 1 m^^s^^, respectively, for temperatures of 1, 2, 
3 and 10 MeV {Bmax can be considered as the limit for infinite temperature). One notes 
that the possibility of bubble nucleation is fairly insensitive to the temperature as soon as 
T exceeds a few MeV, whereas thermally induced bubble nucleation is impossible for T < 2 
MeV (it is known from the stability of ordinary nuclei against decay into quark matter that 
B > (145 MeV)'^). This confirms the estimate in |^. The range of bag constants for which 
a hot neutron star may transform into quark matter is thus roughly 145 MeV < P^/^ < 152 
MeV. 



An interesting feature of the solution is the existence of a maximum in S as a function 
of hb- This indicates that there is a range in densities for which boihng can take place 
for a fixed value of B. If the central density of a neutron star exceeds the upper limit of 
ub permitting boiling, it may therefore happen that boiling is initiated off-center, with po- 
tentially interesting consequences for supernova energetics, neutrino fluxes, gamma-bursters 
etc. This may be explained as follows: With increasing density one has an increase in fi, 
and since Pud ~ At^ while P„ ~ fi^, higher densities must imply still lover Bmax in order 
to satisfy the condition, AP = 0. That the maximum does not occur at the same density 
for all temperatures is due to the /^-dependence of 7. This effect, however, occurs at much 
higher hb and B for the more realistic equation of state discussed below, so we shall not 
pursue the issue further. 

IV. MEAN FIELD APPROXIMATION 

In the following, a mean field model is used to describe the equation of state in the 
hadron phase. The model includes the light hadron octet and is described in detail in Refs. 



Il| , p!2|jl3| . For the quark phase, the only difference from Sec. |T| is that strange quarks are 
introduced in accordance with Eq. (|l^) below. This leads to additional contributions to Eqs. 
(0) and (I). 

Integrating the Oppenheimer-VolkofF equation gives the structure of the initial neutron 
star: Pressure, chemical potential, and number density of each hadron species as a function 
of radius. 

For a flavor conserving phase transformation, the relative number densities, Tj = rii/nB, 
are given by 



(ru\ 

VsJ 



/2 1 1 2 1 1 0\ 

12 10 12 1 
VO 1 1 1 1 2 2/ 






(11) 



The absolute number densities are then given from equahty between the chemical po- 
tential per baryon in the two phases. 

Doing a calculation similar to the one already done for a pure neutron star, one obtains 
a qualitatively similar result with the limits being independent of the nucleation rate con- 
sidered for T > 10 MeV but with slightly higher values of Bmax, and with the maximum in 
B at higher densities. (Fig. ||). 

Again, the temperature variation of the 5-limits show that there is some critical temper- 
ature, Tcrit ~ 3 MeV, below which the conversion takes place only for extremely low values 
of the bag constant. This corresponds to the limit found by Horvath et al. ^ although the 
method for obtaining it is quite different. As mentioned in the previous section, increasing 
the value of T brings the limits closer to the AP = - curve. 

A typical temperature for a newborn neutron star is about 10 MeV at the center with 
off-center temperatures up to 20-30 MeV, but since we have just seen that the exact value is 
unimportant in the high T limit, it is enough to consider only T = 10 MeV (this agrees well 
with the choice of r = 1 s, since a typical cooling time for newborn neutron stars is of this 
order [Q)- The constraints on a 1.4 Mq neutron star is shown on Fig. ^. For B^^'^ < 165 
MeV the center of the star is converted, and for lower values of the bag constant still larger 
fractions of the star will undergo the phase transition in the initial stage of its existence. At 
B^^^ = 145 MeV the entire star is transformed but for this value ordinary nuclei (e.g. ^^Fe) 



would probably be unstable |TH]. For different choices of the hyperon coupling constants. 



the limits on Bmax at the center vary from (162 MeV)^ to (170 MeV)^, with high B^ax 
corresponding to low values of the hyperon-to-hadron coupling constant. Thus a weaker 
coupling of the hyperons tends to destabilize hadronic matter. (The coupling constants 
mentioned above are chosen to fit the A binding energy, as described in Ref. 0]). At larger 
radii, the relative density of strange baryons decreases and the limits converge toward the 
ones shown in Fig. ^. 

So far, we have only considered M = IAMq since observational data seem to suggest 
this as the most typical value. The effect of varying the mass is displayed on Fig. ^. For low 
masses, the star consists almost only of neutrons, even at the center, and thus the deviation 
from (145 MeV)'' seen here is due only to a nonzero pressure, and it appears that in this 
case, a conversion via thermal bubble nucleation is less likely to take place. At masses near 
the maximum mass B^^^ — » 200 MeV so here all but unrealistically high B gives a transition 
into quark matter. 

Turning again to the case M = IAMq, it could be interesting to examine the effects of 
having a non-zero strong coupling constant, a^. For massless quarks, this gives corrections 
to the number density and pressure which to first order is given by Ug = ng^o{l — ^^) and 
Pg = -Pg,o(l — ^)y where denotes the values for a^ = 0. 

Unfortunately, the corresponding expression for 7(as) (and in the case of massive quarks 
also cr{as)) is presently unknown and thus we cannot correctly estimate the effect on the 
nucleation rate. What we can do, however, is to examine the AP = 0-curve, and thereby 
obtain limits on TZ, since this is bound to be below the curve of equal pressure. 

It is seen that although both Bcrit and the limiting bag constant for ud quark matter 
stability {Bcrit taken at r = R corresponding to zero external pressure |]16|) are decreasing 
functions of a^, one has a narrowing of the relevant interval in B (Fig. ^), and above a^ ~ 0.6 
conversion of neutron stars seems to be ruled out. (Similar results were obtained in Ref. 

E3)- 

The effect on Tcru can only be guessed at, but assuming that Wc ~ (1 — ^)", where 
\a\ < 1 — 2, and a^ < 0.6, Tcru should be correct within a factor of 2, and thus the 
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temperatures accompanying supernova explosions should still be enough to ensure conversion 
into strange stars (or hybrid stars) provided that the bag constant is below Bcru- 

V. DISCUSSION AND CONCLUSION 

What we have seen is that if the bag constant lies in the interval where three flavor but 
not two flavor quark matter is stable at zero pressure and temperature (145 MeV < B^^'^ < 
163 MeV, see Ref. |jl5|) then all or parts of a neutron star will be converted into strange 
matter during the flrst seconds of its existence (but note the cautionary remark in [Q). 
The rest will then be transformed either by a slow burning on a time scale of a few seconds 
to a few minutes or by a detonation 0. For bag constants above the stability interval, 
we have seen that a partial transformation is still possible, but since this seems to depend 
heavily on the exact equation of state, one should be careful before drawing any definite 
conclusions. 

Since a large fraction of the star is converted on a relatively short time scale, the released 
energy may well provide a significant contribution to the total energy of a supernova (cf. 
Ref. i). 

Another investigation by Krivoruchenko and Martemyanov Jl^, taking AP = as a 
criterion for a possible transformation into strange stars by a flavor conserving phase transi- 
tion, have found similar results. This is an effect not as much of equivalence of the methods 
used, but rather of the fact, that the high temperatures of newborn neutron stars together 
with the exponential in Eq. (^ causes the rate to be insensitive to T. 

Another interesting feature is that if a star is born with a mass that for a given bag 
constant is too small for the conversion to take place even in the center, then accretion from 
a neighboring star, leading to a higher mass for the neutron star and thus in principle a 
larger transition probability, will only lead to a phase transition via thermal nucleation if 
at the same time the neutron matter is heated to at least 2-5 MeV by the energy released 
by the accretion process or by other mechanisms, such as capture of high energy neutrinos. 



(It is very unlikely that a significant mass can be transferred during the first second or so). 
Thus, one may conclude that if no mechanism for a significant heating of the star can be 
found, the initial mass uniquely determines the future of the star if one has to rely solely 
on thermal nucleation. (Other possible mechanisms that may lead to a transformation were 
mentioned in Sec. |). 

By introducing a non-zero a^, a narrowing in the interesting range for B was seen; both 
as an absolute measure and in terms of the fraction of the interval where uds quark matter 
is stable at zero pressure. 

In this work we have ignored the effect of the mass of the strange quark, which corre- 
sponds to a somewhat inadequate treatment of the equation of state, surface and curvature 
effects. However, even in the center of the star no more than 2-4 % of the quarks in 
the hadrons are s quarks, and thus the effects during bubble nucleation are very small. 
As far as the surface energy is concerned, a more important effect comes from taking 
o'hadron — (SOMcV)'^ (from typical nuclear mass formulae), but even here it turns out that 
47rr^cr ^ 87rrc7, so that inclusion of such a term would not change our conclusions. 
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FIGURES 
FIG. 1. The upper limits on the bag constant allowing boiling of a nonrelativistic neutron gas 

into ud quark matter as a function of the baryon number density in the hadron phase is shown for 

different nucleation rates and temperatures. The upper curve for each temperature corresponds to 

a rate of one nucleation per km per Gyr; the lower to one per m^ per second. As comparison is 

shown the AP = line for a Bethe- Johnson equation of state (BJ). This is seen to deviate only 

for very large densities. 

FIG. 2. Limits on the bag constant in the mean field approximation. Notation as in Fig. |^. 
(The model used correspond to x=0.6 in Ref. p3[|). 



FIG. 3. Limits on the bag constant for a 1.4 Mq neutron star. Notation as in Fig. ||. 

FIG. 4. Critical bag constants in the center of neutron stars as a function of stellar mass. 
Notation as in Fig. |l[ 

FIG. 5. Curves representing pressure equilibrium between the two phases for various values of 
the strong coupling constant, a^. 
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